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Abstract 

We derive model-independent, universal upper bounds on the Operator Product Expan- 
sion (OPE) coefficients in unitary 4-dimensional Conformal Field Theories. The method 
uses the conformal block decomposition and the crossing symmetry constraint of the 4-point 
function. In particular, the OPE coefficient of three identical dimension d scalar primaries 
is found to be bounded by ~ 10(d — 1) for 1 < d < 1.7. This puts strong limits on unparticle 
self-interaction cross sections at the LHC. 
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In this paper we will answer, in a particular well-defined context, the question: Is there an 
upper bound to the interaction strength in relativistic Quantum Field Theory (rQFT)? 

Intuitive reasons suggest that such a bound exists. Take QCD as a representative real-world 
example. At energies E above the scale Aqcd ~ 1 GeV, this is a perturbative theory of interacting 
quarks and gluons, and the interaction strength is measured by the dimensionless running coupling 
g s {E). The coupling starts small at very high energies E 3> Aqcd and grows at low energies, 
formally becoming infinite at E ~ Aqcd- However, perturbative expansion breaks down before 
this happens. L-loop diagrams are suppressed by factors ~ (g 2 /167r 2 ) i . As soon as g s ~ 4tt, 
all loop orders contribute equally. Thus in perturbation theory it is impossible to get couplings 
stronger than about An. 

To recall what happens beyond perturbation theory, let us look at the same theory at energies 
below Aqcd- In this regime the appropriate degrees of freedom are hadrons, and their interactions 
can be described by an effective lagrangian. For instance, pion-pion scattering at low energies is 
described by the chiral lagrangian 

C = ^Tr|^f/| 2 + ..., U = ex P (t2n a T a /f 7T ), 

where f n ~ 93 MeV is the pion decay constant, T a are the SU(2) generators and . . . stand for the 
chiral symmetry breaking terms. The dimensionless quartic pion coupling can be defined from the 
2 — > 2 scattering amplitude; it grows with energy as A ~ (E/ f n ) 2 . If the chiral lagrangian is valid up 
to energies ~ Aqcd and is stable under radiative corrections, we should have A(Aqcd)/16vt 2 < 1, 
or Aqcd ^ Anf^. Experimentally this bound is satisfied and, moreover, near-saturated. This 
observation forms the basis of the Naive Dimensional Analysis [1] method of estimating couplings 
in strongly coupled theories. 

While the above arguments are appealing, at present it is unknown if they can be turned into 
a theorem, or even how to formulate such a general theorem. In order to make progress, in what 
follows we will assume that we have a Conformal Field Theory (CFT), i.e. an rQFT invariant 
under the action of the conformal group [2]. 

CFTs form an important subclass of rQFTs. Presumably, any unitary, scale invariant rQFT 
is conformally invariant. This is proved in D = 2 spacetime dimensions under very mild tech- 
nical assumptions [3J,[4], and no counterexamples are known in D > 3. Since scale invariance 
is ubiquitous (think of any RG-flow fixed point), this would make conformal invariance equally 
ubiquitous. Unitarity is however crucial here: without unitarity simple physical counterexamples 
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exist, e.g. theory of elasticity [BJ. We are interested in applications to particle physics, thus we 
will assume unitarity, and will work in D = 4. 

There are many known or conjectured classes of four-dimensional CFTs. For example, M = 1 
supersymmetric QCD with N c colors and Nf flavors flow to a CFT in the infrared as long as 3/2 < 
Nf/N c < 3 [7J. Large N c analysis [8] and lattice simulations [9] suggest that a similar 'conformal 
window' exists also without supersymmetry. Another famous example is the TV = 4 super Yang- 
Mills (SYM), which is conformal for any coupling and any N c . At large 't Hooft coupling and 
large N c it can be described via the AdS/CFT correspondence. Many deformations preserving 
conformal symmetry are known on both field theory and gravity sides of the correspondence [11] . 
Our discussion will be general and will in principle apply to all the above examples. 

The D = 4 conformal group is finite dimensional; it is obtained from the Poincare group by 
adding the generators of dilatation T> and of special conformal transformations K,^. The local 
quantum fields 0(x) are eigenstates of V, [V, O(0)] = zAO(O), where the eigenvalue A is called 
the scaling dimension. The /C M acts as a lowering operator for the scaling dimension, and the 
corresponding 'lowest- weight states', fields satisfying [/C M , O(0)] = 0, play a special role. They are 
called primaries. All other fields can be obtained from primaries by taking derivatives and are 
called descendants. 

Conformal symmetry constraints the 2- and 3-point functions of primary fields to have partic- 
ularly simple form. For scalar primaries, we have: 



Eq. (II]) says that a diagonal basis can be chosen in the space of primary fields, and sets the 
normalization. Eq. ^ then defines coefficients c^fc. These same coefficients appear in the Operator 
Product Expansion (OPE) 



where . . . stands for the contributions of higher spin primaries and of descendants. 

In CFT, any n-point function can be, in principle, reduced to a sum of products of 2-point 
functions by repeated application of the OPE, with coefficients given by products of c^'s and 
of their higher spin generalizations. In this sense, cafe's play in CFT a role similar to that of 



(O i (x 1 )O j (x 2 )) = 6 ij (x 2 12 )- A , 

(Oi{x x )o s {x 2 )o k {x*)) = c ijk (xl 2 r^(x 2 13 y^(x 2 23 y^ , 

x\ = (xi - xj) 2 , p ijk = (Aj - Aj - A fe )/2 . 



(1) 
(2) 



0,(X)0,(0) ~ (^)-(A I+ A,)/2 ^ + c .. k tf)W 0h QS) +...}, 
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the (renormalized) coupling constants in perturbation theory, measuring interaction strength. We 
thus have the following CFT version of our initial question: Is there an upper bound to the OPE 
coefficients, valid in an arbitrary unitary CFT in D = 4 ? We will now proceed to show that such 
a universal bound indeed exists. 

Let us pick a hermitean scalar primary of scaling dimension d and consider its OPE with 
itself: 

0(*)0(O)~(*V{l+ J2 E c AlJ (x 2 )( A -^V 1 ---^O w ... w (0) +...}. (3) 

1=0,2,4... A>A min (0 

This time we show explicitly contributions of both scalars (I = 0) and of higher spin primaries 
Op-L^m which are symmetric traceless tensors. Spin I has to be even by the Bose symmetry. Lower 
bounds on the dimension A of a spin I primary: 

A min (/ = 0) = 1, A min (Z > 1) = I + 2 , 

(unitarity bounds) are known to follow from unitarity [12J. Only special fields may saturate these 
bounds: a free scalar (1 = 0), conserved currents (1 = 1), and the stress tensor (1 = 2). Higher I 
conserved currents, present in free theories, also saturate the bounds. 

An interesting object to study is the 4-point function of <fi, constrained by conformal symmetry 
to have the form 

q(u, v) 

(<f)(x 1 )<j)(x2)<f)(x 3 )(f)(x4 : )) = ' 2d , (4) 

x 12 x u 

where u = x\ 2 x\J (x\ z x\^) , v = x\^x\^j '(xf 3 a;| 4 ) are the conformal cross-ratios. The same 4-point 
function can be reduced to a sum of 2-point functions by applying the OPE in the 12 and 34 
channels. Cross terms of different primary families drop out of this representation because of 
Eq. and its higher spin analogue. Terms involving the same primary and its descendants can 
be resummed in closed form. As a result, we get the conformal block decomposition 

g(u,v) = 1 + ^2pa,i9a,i(u,v) , Pa,i = c ai , (5) 

where [13] 

(-V zz 

9a,i(u,v) = — z [k A+ i(z)k A -i- 2 (z) -(z++ z)\ , 

l l z — z 

k p (x)= x^SF^P /2,(3/2, (3;x) , 
u = zz, v = (1 — z)(l — z) . 
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This decomposition is expected to converge at least in the circle \z\ < l,\z\ < 1, which corresponds 
to being able to fit a sphere centered at x\ which separates 22 from 23 and £4 [14J. 

The 4-point function Q must be crossing-symmetric under the x\ <-> X2 and x\ <-> X3 ex- 
changes. The first crossing is manifest since only even spins contribute to the OPE. The second 
one gives a nontrivial constraint 

v d g(u, v) = u d g(v, u) . (6) 

Decomposition (J5]) must be consistent with this constraint. Separating the contribution of the 
unit operator, we obtain the sum rule 

1 = 22PA,iFd,A,i(u,v) , (7) 
j? , \ v d g Aj i(u,v) -u d g Ad (v,u) 

F d ,AA u i v ) = a d • 8 

u a — v d 

As we will now show, this equation can be used to get an upper bound on c A ,i- 

Crucially, coefficients c A ,i are real, and thus p A j > 0. This can be related to the absence 
of parity violation in the conformal 3-point function of two scalars and a symmetric tensor [15] . 
Eq. ([7| then allows a geometric interpretation: when p A i > are allowed to vary, the RHS fills 
a convex cone Cd in the vector space V whose elements are two-variable functions. We say that 
this cone is generated by functions Fd,A,i{u, v), I = 0, 2, 4, . . . , A > A min (7). Eq. ^ expresses the 
fact that the function f(u,v) = 1 belongs to this cone. 

Let us now pick a particular field O^j and rewrite (l7|) as 



1 - PA,iF d ,A,i{u, v) = PA,iF dA ,i (u,v). (9) 

As J is increased, the vector corresponding to 1 — Pa jF d A i( u i v ) moves in the vector space. 
Suppose that for all p A j above some critical value p cr this vector stays out of the cone Cd- Then 
p cr provides a bound on the squared OPE coefficient |c A /| 2 . This bound will depend on d,A,l, 
but will be valid in any unitary CFT. 

To find p cr , we employ the method of linear functionals developed in [T5].[16j. Recall that a 
linear functional is a linear map A from V to real numbers: 

A : V - R, AfoiFi] = OiA[Fi] . (10) 

Suppose that we found a functional which is positive on all functions generating the cone Cd'- 

A[F, A/ ]>0. (11) 
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We will normalize this functional by the condition 



A[l] = 1. (12) 
Since for such A Eq. ([7]) implies A[l — PajF^aj] > 0> we would get an upper bound: 

PA,l<Pcr(A) = l/A{F dAI \. (13) 



To make this bound as strong as possible, we will impose, in addition to (11 ), (12), an extremality 
condition 

MF d ,A,i\ - m ax . (14) 

We will use linear functionals given by a finite linear combination of derivatives evaluated at 
a given point. More precisely, we will use functionals of the form: 

A[F]= KmF i2n > 2m \ N = 3, (15) 

n,m> 0, n+m< N 
F (2n,2m) = tfntfm 

^ = l/2 + a + 6, z = 1/2 + a- 6. (17) 

Here A n>m are fixed real numbers defining the functional. The symmetric point a = b = is chosen 
as in [I5],[T6] since the sum rule is expected to converge fastest here, and because the functions 
Fd,A,i are even in both variables with respect to this point. This is why only even-order derivatives 



are included in (15). 



Eqs. (11), (12), (14) define an optimization problem for the coefficients A nm . The constraints 



are given by linear equations and inequalities, and the cost function is also linear, which makes it a 



linear programming problem. Although the number of constraints in (11) is formally infinite, they 
can be reduced to a finite number by discretizing A and truncating at large A and /, where the 
constraints approach a calculable asymptotic form. The reduced problem can be efficiently solved 
by well-known numerical methods, such as the simplex method. A found solution can be then 
checked to see if it also solves the full problem. This procedure was developed and successfully 
used in a different context in H51.I16I. 



In this work, we used this procedure to compute bounds on the OPE coefficients c^o when O 
is a scalar field (1 = 0). We will now present our numerical results. Fig. 1 concerns the case when 
the dimension of <fi is close to that of a free field, 1 < d < 1.1. Notice the bell-shaped form of the 
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max |c #(9 | 




Figure 1: Theoretical upper bound for the OPE coefficient c^o a function of the 
dimension A of the scalar field O. The curves correspond to the <ft's dimension fixed 
at d = 1.005, 1.02, 1.05, 1.1 (from below up). The bound was computed for each of the 
shown points, and the curves in between were obtained by interpolation. 

bound, peaked at A ~ For d — > 1 the bound evidently tends to zero everywhere except near 
A = 2. This means that the free field theory limit is approached continuously: for d = 1 the only 
scalar operator in the 0x0 OPE is the : 2 : of dimension 2. In Fig. 2 we present a similar plot for 
1.2 < d < 1.7. Notice that the bounds in Figs. 1,2 go to zero as A — > 1. This is expected in view 
of the general theorem that a dimension 1 scalar must be free, hence decoupled from everything 
else in the CFT. 

A text file with the coefficients of the linear functionals used to obtained the bounds plotted 
in Figs. 1,2 is included in the source file of this arXiv submission. The reader may check that they 



indeed satisfy the constraints (11). 



We have only explored the range d < 1.7 for the following reason: starting from d ~ 1.75, we 



found that there is no functional of the form (15) satisfying the constraints (11), (12). We expect 



that a bound exists also for larger d, but to find it one needs to use more general functionals, 
e.g. involving more derivatives (i.e. with higher N). This will also give improved bounds in the 

lr This shape makes it tempting to draw an analogy with the Breit-Wigner formula, especially since the dilatation 
operator T> plays the role of energy in radial quantization. 
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max |c #c | 




A 

Figure 2: Same as Fig. 1 for the 's dimension fixed at d = 1.2 , 1.3 , 1.4 , 1.5 , 1.6 , 1.7 
(from below up). 



max \c m \ 




Figure 3: Theoretical upper bound for the OPE coefficient c^ as a function of <j)'s 
dimension d. 
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range of d that we considered. This is left for future work. 

On the other hand, the restriction to 1 < A < 3 in Figs. 1,2 is not essential: our method would 
also give bounds beyond this range. In fact, any of the functionals derived for 1 < A < 3 could 
be used to compute a sub-optimal but valid bound for larger A (as well as for I > 0) via Eq. (13). 

It would be interesting to study the asymptotic behavior of the bound at large A. A conser- 
vative upper estimate can be obtained from the known asymptotics of F d ^j and its derivatives 
|15j . if we assume that the functional A in (13) is A-independent. This way one concludes that 
the bound cannot grow faster than exponentially: {c^ol = 0(q A ), q = (v^2 + l)/2. However, this 
is likely an overestimate, since the optimal functional A, as determined by Eq. (14), will likely 
depend on A. 

It would be also interesting to derive analogous bounds in two spacetime dimensions, where 
explicit expressions for conformal blocks are also known [13J. 

As a phenomenoligical application of our results, consider the unparticle physics scenario [17] . 
Unparticle self-interactions were considered in [18] (see also [19]) a prominent feature of such 
scenarios, giving rise to processes like gg — > — > <f><fi — > 47. The cross section for this process is 
proportional to the square of the self-coupling OPE coefficient c^, where is a scalar operator 
from a hidden-sector CFT (unparticle) with non-renormalizible couplings to gluons and photons. 
In [18] , the values of these coefficients were kept as arbitrary parameters, unconstrained by prime 
principles, and only experimental constraints from the Tevatron were imposed, which led to a 
possibility of spectacularly large cross sections at the LHC. In Fig. 3 we plot our theoretical 
upper bound on (extracted from Figs. 1,2 by setting A = d). The values of used in [IB] 
exceed our bound by 2 4 orders of magnitude]^] We conclude that a revision of the studies in 
[18], [20], taking into account our bounds, is necessary. 

As a purely field-theoretical application, consider the Af = 4 SYM theory already mentioned 
above, which is conformal for any value of the 't Hooft coupling A = gy M N c . The region of small A 
is accessible via perturbation theory, while large A (and large N c ) are accessible via the AdS/CFT 
correspondence. Moreover, the large N c theory is integrable, which allows to interpolate between 
the two regimes and perform various nontrivial checks [21]. As A is increased from to 00, the 
spectrum of the theory changes, and anomalous dimensions of some local fields are certain to 
become large. For example, at large N c the fields which do not map onto supergravity modes 



2 For proper comparison note that the normalization of the unparticle OPE coefficients C'd used in |18j is related 
to our normalization via Cd — g^ 2 (\^<i\/9d) 3 ^ 2c 4>4>4> where Bd is given in [TS] and gd = 4 2_d 7r 2 r(2 — d)/T(d). 
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on AdS$ x S 5 have anomalous dimensions growing for large A as A 1/4 [22]- Now one could ask 
what happens to the OPE coefficients, whether they can have similar growth. From our results, 
assuming that they can be extended to d > 1.7 as discussed above, it follows that no matter how 
large A is, the OPE coefficients of fields with low dimensions will stay bounded. It should be 
noted that this conclusion is nontrivial only for small N c , since at large N c the OPE 0\ x O2 is 
known to factorize, with the composite "multi-trace" fields : O1O2 '■ appearing with the coefficient 
1 + e>(l/iV c 2 ) while all other fields 1/N C suppressed [23]. 

In summary, we have presented theoretical upper bounds on the OPE coefficients of two 
identical scalars and a third scalar, valid in an arbitrary unitary CFT. Our results are based on 
imposing crossing symmetry on the conformal block decomposition of a scalar 4-point function. 
They imply that, in a certain sense, interaction strength remains limited even in theories like 
M = 4 SYM (or its many known conformal deformations) where a coupling A can be taken to 
infinity. They also lead to strong bounds on the cross sections of unparticle self-interaction-type 
processes at future colliders. 
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